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ABSTRACT 

We perform a duality transformation that allows one to express the partition func- 
tion of the ci-dimensional Ising model with random nearest neighbor coupling in terms 
of new spin variables defined on the square plaquettes of the lattice. The dual model 
is solved in the mean field approximation. 
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The Ising model with random coupling plays a fundamental role in the theory of 
disordered systems. In this field, one of the major results is the Parisi solution of the 
infinite range model where at low temperature the system becomes a spin glass with a 
replica symmetry breaking [1]. However, there is no exact solution for the model with 
nearest neighbor interaction, and it is still unclear whereas a glassy phase is present 
in three dimensions at finite temperature. 

In this letter, we perform a duality transformation of the Ising model with random 
nearest neighbor coupling that assume the values = ±1 with equal probability. The 
model is thus defined on a dual lattice where the spin variables are attached to the 
square plaquettes. The advantage is that the non-linear part of the dual Hamiltonian 
has constant coefficients instead of random ones. It is therefore possible to use the 
standard mean field approximation to estimate the quenched free energy. However, in 
the dual lattice the ratio between number of spins and number of links increases with 
the dimension at difference with what happens in the original lattice so that the mean 
field approximation becomes worst at increasing the dimensionality Our solution is 
thus optimal in two dimensions where it gives a rather good estimate of the ground 
state energy. 

The partition function of the d-dimensional Ising models on a lattice of N sites 
with nearest neighbour couplings Jjj which are independent identically distributed 
random variables, in absence of external magnetic field, is 

Zn(0, {Jij}) = II exptfJijWj) (!) 

where the sum runs over the spin configurations {a}, and the pruduct over the nearest 
neighbor sites One is interested in computing the quenched free energy 



where A indicates the average of an observable A over the distribution of the random 
coupling. The quenched free energy is a self-averaging quantity, i.e. it is obtained 
in the thermodynamic limit for almost all realizations of disorder [1]. Even in one 
dimension, it is difficult to find an exact solution for / in presence of an external 
constant magnetic field [2]. On the other hand, it is trivial to compute the so-called 
annealed free energy 

corresponding to the free energy of a system where the random coupling arc not 
quenched but can thermalize with a relaxation time comparable to that one of the 
spin variables. An easy calculation shows that in our case 

fa = -i(ln2 + dlncosh/3) (4) 

However, f a is a very poor approximation of the quenched free energy, and is not able 
to capture the qualitative features of the model. 

In order to estimate (1), it is convenient to use the link variable Xy = Oi <jj, since 
only terms corresponding to products of the variables x^ on close loops survive after 
summing over the spin configurations: on every close loop of the lattice Y\xij = 1, 
while Y\ x ij = Ca o& for a path from the site a to the site b. A moment of reflection 
shows that it is sufficient to fix Y\xij = 1 on the elementary square plaquettes V to 
automatically fix it on all the close loops. The partition function thus becomes 

W{J«})=E II V 1 U iJ ( 5 ) 

where the number of plaquettes is N p = d(d — l)N/2, and we have introduced the 
plaquette variable Xi = Y\ v x^. 

For dichotomic random coupling Jjj = ±1 with equal probability, the free energy 
of the model is invariant under the gauge transformation x^ — > Jjj Zy , so that one 
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has 

^ = e n n ^ ( g ) 

where Jj = J|-p J»j i s again a dichotomic random variable (the 'frustration' [3] of the 
plaquette Pi). It is worth remarking that (6) gives the partition function in terms of a 
sum over the 2 dN configurations of the independent random variables = ±1 with 
probability 

In the following we shall indicate the average of an observable A over such a normalized 
weight by (A), e.g (x^) = tanh/3 and (a?j) = tanh 4 /3. With such a notation, the 
partition function assumes the compact form 

N p 

Z N = 2<< W - JV *> cosh dA, (/3) ( J](l + xJi) ) (8) 

i=l 

Now comes the key step. We estimate the average in (8) by a geometrical construction. 
Let us introduce the dual lattice [4] as the lattice whose sites are located at the centers 
of each square of the original lattice. A dual spin variable is attached to each square 
plaquette and can assume only the values <r, = ±1 with equal probability, so that one 
has the identity 

(l + xj i )= E {Xi J,) (1+ ^ )/2 - (9) 

Since there is a one-to-one correspondence between links on the original and on the 
dual lattice, we can compute the link-average noting that 

<n-f^> = ( n *F w(4j) fid?fci > = n ^-^m m 

where Efce(i j) i s the sum on the 2(d — 1) nearest neighbor plaquettes variables with 
a common link (i, j) and for sake of simplifying the notation we have defined the new 
link variable 

p$' j) = n (") 

ke(i,j) 
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Inserting (10) and (9) into (8) one has 



Z N = 2^ dN - N ^ cosh dJV (/?) LI ^ (1+ffs)/2 II ePV*-V (12) 

{a} i=1 

where the inverse temperature of the dual model is 

= - ^ In tanh(/3) (13) 

that vanishes when the temperature T — /3 _1 — > 0. The quenched free energy (2) thus 
is 

-/W) = f ((2-d)ln2 + lnsinh(2/3)-(rf-l)^/(/3)) (14) 
where the free energy of the dual model is 

70) = - lim ^-ln^ p (15) 



with 



^ P = E ex p ^E ^ <J) I (-i)^>' {7i}) (i6) 

{<r,} \ W) 



and 



1 + a% 



Since the Hamiltonian of the dual model is defined by the relation Z^ p = ^{^} e 
from (16) one sees that 

tf = -E p <P - E ln ( J «) ( 18 ) 

Let us stress that the non-linear term of the dual Hamiltonian (18) does not depend 
on the random coupling like in the original model and the randomness enters via the 
presence of a sign. In fact, the weight exp(—/3H) does not defines a standard Gibbs 
probability measure on the dual lattice: it defines a signed probability measure. For 

(i j) ~ ~ 

instance, in 2d, P d = o% <Jj so that the Gibbs measure of a configuration of the 
dual random Ising model differs from that one of the pure Ising model only for the 
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presence of the random sign determined by the index K related to the frustrations of 
the square plaquettes { Ji}. 

This is the first result of our letter. Its importance stems from the fact that it is 
now possible to linearize the Hamiltonian (18) by introducing the magnetization 



m= lim — y^CTi (19) 

y i 

Indeed, if we neglect the fluctuations, the non-linear term of (18) can be estimated as 

so that (16) becomes 

Z Np =J2 n^ (1+ ^ )/2 exp(^d7Vm 2 ( d - 1 )) (20) 

Let us find the mean field solution, by using an auxiliary field Recalling the saddle 
point method one immediately sees that in the limit N — > oo, 



,Nd{} m >i<>-» „ d<I>e X p(Npd(C d m<S> 2d - 3 - $ 2 ( d " 1 ))) (21) 

J — oo 



where the constant d is determined by maximizing the argument of the exponential 
and reads 

C d = 2 (d - 1) ( 2d - 3 )-(d-3/2)/(d-i) (22) 
As a consequence one can write the partition function as 

Z Np ~ max exp (-A^ 2 ^" 1 ) ) £ [] jf + ^ }/ V^* 2 "^ (23) 

Now the effective Hamiltonian is linear and one can explicitly carry out the sum over 
the N p independent random variables <7j, obtaining 

-0f0) =dmax| ^-—^ ln(2sinh( 7d /3$ 2<i - 3 )) - /jf* 2 ^ -1 ^ (24) 
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with -fd = 2 Cd d. The maximum is realized by the value of $* that is the solution of 
the self-consistency equation 

{2d - g 3hd coth( 7 ^$ 2d - 3 ) = $ (25) 

In 2d, (25) assumes the simple form coth(8/3<I>) = <f>. The graphical solution of this 
implicit equation is showed in fig 1. One sees that should be always larger than 
unity and at — > oo (infinite temperature T = /3 _1 limit) <&* = 1. It can appear 
rather odd that in the dual model the magnetization $* > 1. This stems from the 
fact that the Gibbs probability measure exp(—0H) is a signed measure because the 
random coupling is transformed into a complex random magnetic field in (15). From 
fig 1, it is also clear that the mean field solution does not exhibit phase transitions at 
finite temperature. However there is an essential singularity at T = 0, since inserting 
(24) into (15) and (14) one sees that / - exp(l/T) for T -> 0. 

It is important to stress that the mean field solution does not improve at increasing 
the dimensionality, since the ratio N p /N between number of plaquette spins and 
number of links in the dual lattice increases as (d — 1) at difference with the standard 
Ising model where the ratio of spins over links decreases as d^ 1 . 

It is possible to explicitly solve the self-consistency equation for — > where (25) 
becomes 

$* - (12/3)~ 1/2 ford = 2 (26a) 



M _ 3X l/(2 d - 2 ) 



and 

$* „ i^^J £-V(2d-2) for d > 3 (266) 

Such a relation shows that when d — > oo one has $* — > 1 for T — > and then for 
all T's. The high dimension limit is therefore trivial. The mean field approximation 
works at its best in two dimensions. For instance, the zero temperature energy of the 
mean field solution is E a = — 3d/ A which is a fair estimate in 2d where the numerical 
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simulations [5] give E = —1.404 ± 0.002. In fig 2, we show the free energy as a 
function of T in 2d compared with the annealed free energy (4). One sees that entropy 
is negative at low temperature, thus indicating that the solution is unphysical. As a 
consequence, a better estimate of the ground state energy is given by the maximum 
of /(/3), following a standard argument of Toulouse and Vannimenus [6], and one has 
E a < max f(P) = -1.468. 

In conclusion, we have obtained two main results. (1) Formulation of the random 
coupling Ising model on the dual lattice made of square plaquettes. The dual model has 
signed Gibbs probability measure and magnetization larger than unity. (2) Application 
of the mean field approximation to solve the dual model. The approximation is sensible 
at low dimension. 

In our opinion, the mean field approach is very promising at least in two 
dimensions. It has a good heuristic power and there are still a lot of open problems in 
its framework, such as finding a Ginzburg-Landau criterion or refining the mean field 
approximation in a cluster expansion scheme. 

We acknowledge the financial support of the I.N.F.N., National Laboratories of 
Gran Sasso. 
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Figure Captions 

1 Graphical solution of the implicit equation (23) in 2d, at T = = 1 
corresponding to (3 = 0.136... The full lines are coth(8/3$) versus $ and the 
straight line $ = 

2 Random Ising model in 2d: the annealed free energy f a (dashed line) and the mean 
field solution (full line) versus temperature T — fi^ 1 . The dashed lines are the 
Maxwell constructions obtained by imposing that the free energy is a monotonous 
non-decreasing function of T. The annealed solution gives a ground state energy 
E = —1.559; the mean field solution gives E = —1.468; the numerical result of 
[5] is E a = -1.404 ±0.002 
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